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Chiral phase structure of QCD with many flavors 
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We investigate QCD with a large number of massless flavors with the aid of renormalization 
group flow equations. We determine the critical number of flavors separating the phases with and 
without chiral symmetry breaking in SU(N C ) gauge theory with many fermion flavors. Our analysis 
includes all possible fermionic interaction channels in the pointlike four-fermion limit. Constraints 
from gauge invariance are resolved explicitly and regulator-scheme dependencies are studied. Our 
findings confirm the existence of an Nf window where the system is asymptotically free in the 
ultraviolet, but remains massless and chirally invariant on all scales, approaching a conformal fixed 
point in the infrared. Our prediction for the critical number of flavors of the zero-temperature chiral 
phase transition in SU(3) is Nf 1 = 10.0 ± 0.29(fermion) ^J'gg (gluon), with the errors arising from 
approximations in the fermionic and gluonic sectors, respectively. 
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I. INTRODUCTION 

In order to investigate the dynamics of strongly in- 
teracting gauge systems, it has often been a successful 
strategy to deform the desired system into a more acces- 
sible one. For instance, the approximation of a contin- 
uum gauge theory by a finite lattice gauge theory can be 
viewed as such a deformation. In the same spirit, the 
addition of more symmetries as in supersymmetric ver- 
sions of gauge theories represents such a more tractable 
deformation. In the present work, we consider the de- 
formation of SU(N C ) gauge theories by many massless 
fermion flavors. Since light fcrmions have the perturba- 
tive tendency to screen long-range forces, a large fermion 
flavor number Nf has the potential to move the system 
towards weaker gauge interaction for which analytic tools 
are more powerful. 

The relevance of fermionic screening at large Nf is ul- 
timately observable for Nf > Nf" ' :— ^N c where asymp- 
totic freedom is lost. But even for smaller Nf, fermionic 
screening is first signaled by the second (3 function co- 
efficient of the gauge coupling, reversing its sign for 



N f > 



34Ng 
13N*- 



This sign change induces a second zero 



of the perturbative j3 function, implying an infrared (IR) 
attractive fixed point of the gauge coupling a* > 0. As 
was first argued in [| , the fixed-point value a* is small for 
Nf < Nf ' f ', supporting the possibility of a perturbative 
analysis. If this fixed point persists on all scales, the sys- 
tem approaches a conformally invariant limit in the deep 
IR, keeping massless quark and gluon excitations in the 
spectrum. This scenario can collapse owing to nonpertur- 
bative phenomena such as the spontaneous break-down 
of chiral symmetry. The latter renders the fermions mas- 
sive, implying their decoupling at low scales. Fermionic 
screening properties are thus switched off in the deep IR, 
and the system is characterized by strongly coupled glue 
and the Goldstone bosons of chiral symmetry breaking 
(XSB). 

Since xSB is triggered by the strength of the gauge 



interactions which in turn depends on the amount of 
fermionic screening, we expect the conformal scenario to 
hold for large Nf sufficiently close to Nf' '. The broken 
phase is supposed to occur for smaller Nf with a critical 
Nf r separating the two phases. In fact, evidence for this 
phase structure has been collected from various methods, 
including (improved) ladder truncated Dyson-Schwinger 
equations 0, .lattice QCD [jj, anomaly-induced ef- 
fective potentials and instanton models Q, with es- 
timates for the critical number of flavors ranging from 
N^ r ~ 5 to N^ r ~ 12 for SU(3) gauge theory. 

An analysis of the quark-scattering amplitude using 
the functional methods 0, particularly reveals that 
the nature of the phase transition, though continuous, is 
not conventionally second order. 1 This is most promi- 
nently visible in the fact that there appear to be no light 
scalar states in terms of which an effective theory could 
be constructed on the symmetric side (Nf > Nf r ) of the 
phase transition 0, 0, EJ Any attempt at construct- 
ing a Ginzburg-Landau-type of effective potential for the 
chiral order parameter hence produces discontinuities in 
the potential parameters. On the other hand, the or- 
der parameter in the form of the chiral condensate, i.e., 
the minimum of this would-be potential, changes contin- 
uously across the phase transition. 

In this work, we analyze the phase structure of many- 
flavor QCD with the aid of the functional renormal- 
ization group (RG). Beyond its conceptual simplicity, 
the approach operates in continuous spacetime, supports 
an explicit implementation of chiral symmetry, and in- 
cludes a rcsummation of beyond-ladder diagrams already 
in simple truncations. Moreover, we monitor gauge in- 
variance with the aid of Ward-Takahashi identities and 
study regulator-scheme dependencies. Since the critical 



1 Owing to the different conformal properties of the system on 
each side, the transition is often referred to as a conformal phase 
transition 0. 



2 



number of flavors is a universal quantity, the scheme- 
dependencies are a direct measure of the approximation 
uncertainties of our nonpcrturbative truncation. Apart 
from a quantitative prediction for N", our results re- 
veal further facets of the nature of the phase transi- 
tion. Across the phase transition, we particularly observe 
discontinuities in the effective fermionic self-interactions 
which are generated by the RG flow. 

For our analysis, we employ the functional RG formu- 
lated in terms of a flow equation for the effective average 
action T k 



d t T k = -STrd t R k (if 5 



RkY 



t = In- 



A 



uv 



(1) 



serving as an alternative definition of quantum field the- 
ory. T k is a free-energy functional that interpolates be- 
tween the bare action Ffc = A uv = S and the full quantum 
effective action T = r k= o- Here, R k denotes a regulator 
function that specifies the details of the momentum-shell 
integrations, the variation of which will provide us with 
an error estimate of our approximations. We devote Sect. 
|Il|to a discussion of our approximation scheme, which 
represents a truncation of the full theory to the relevant 
operators for the present problem. Results and conclu- 
sions are presented in Sects. IIIII and ITVl 



ghost sector is tacitly assumed). The gauge-field kinetic 
term is accompanied by a fc-dependent wave-function 
renormalization Za, the fermionic one by Z^,. Similarly, 
Z\, Z+, Z_, Z G and Zy\ are the vertex renormalizations, 
whereas g, A denote the bare couplings. Renormalized 
dimcnsionlcss couplings arc defined by 



9Z1 



^ A V> 



Zik Xi 



(3) 



We work in the Landau gauge, a = 0, which is known to 
be a fixed point of the renormalization group ^3 an d has 
the additional advantage that the fermionic wave func- 
tion is not renormalized in our truncation, such that we 
can choose Z$ = 1. 

The four-fermion interactions can be classified accord- 
ing to their color and flavor structure. Color and flavor 
singlets are 



(V-A) 
(V+A) 



(?7^) 2 
(V^) 2 



(V'7m75V') 5 



(4) 
(5) 



where color . . . ) and flavor (a, 6, ... ) indices are con- 
tracted pairwise, e.g., (ipi/)) = (ip^f)- The remaining 
operators have non-trivial color or flavor structure 



II. FLOW EQUATIONS 

Our primary aim is a reliable determination of the crit- 
ical flavor number N". For this, the approach of the 
phase transition from the symmetric side is technically 
advantageous, since the relevant degrees of freedom are 
expected to remain the same during the flow from the 
ultraviolet (UV) to the infrared (IR): quarks and gluons. 
Therefore, we solve the flow equation JTJ in a truncated 
subspacc of all possible action functionals. 

In addition to standard gauge and fermion sectors, 
we include all possible pointlike four-fermion interactions 
that are compatible with SU(N C ) gauge symmetry and a 
chiral SU(Nf)L x SU(Nf)R flavor symmetry for Nf fermion 
species 2 , 



r fc = 



WZ^ + Z 1 gM+^Fr F U 
1 
2 
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Z_A_(V-A) + Z+X + (V+A) + Z CT A CT (S-P) 
+Z VA AvA[2(V-A) ad j+ (1/N C )(V-A)]1 . (2) 



(S-P) 
(V-A) adj 



(^ V) 2 - (V/ 7 5^ b ) 2 (6) 
(V7^V) 2 + W'7m7 5 T» 2 , (7) 



where (ip a, 4 )b ) 2 = i?^ 1 ^ ip a , etc., and (T z )ij denotes the 
generators of the gauge group in the fundamental repre- 
sentation. The set of fermionic self-interactions occurring 
in Eq. @ forms a complete basis. Any other pointlike 
four-fermion interaction invariant under SU(N C ) gauge 
symmetry and SU(Nf)L x SU(Nf)R flavor symmetry is 
reducible by means of Ficrz transformations. 

It is important to stress that the effective action © is 
in the QCD universality class 0], as long as the four- 
fermion interactions are small in the UV and thus RG 
irrelevant. Only nonperturbativcly large (^Ai)'s would 
bring us into another universality class with x^B of 
Nambu-Jona-Lasinio (NJL) type, which will not be con- 
sidered in this work. 

Plugging the truncated effective action (J2J) into the 
flow equation, we obtain the follo wing (3 functions for 
the dimcnsionlcss couplings Xi (see 



ving / 
H): 



Here = A Z T Z , F^ v = F Z V T Z denote the nonabelian 

gauge potential and field strength (the inclusion of a =2X — 4v ^ FB )' 4 



N, 



-,g 2 A_ - 3g 2 A V A 



(8) 



2 We note that only the four-fermion interactions arc manifestly in- 
variant under local gauge transformations for all possible choices 
of the couplings. Gauge invariance of the remaining terms is 
governed by (modified) Ward -Takahashi identities as discussed 
for the present system in 1 1 ill . 
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12 + 9N 2 
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-M F 5' 4 {-N f N c (A?. + A5.) + Ai 

-2(N C + Nf)A_ Ay A + NfA+A, + 2A 2 /A }, 
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d t X+= 2A 4 



A ;( FB )> 4 



N, 



(9) 



(FB),4 



12 + 3N 



K, 



3Aj 



2N c N f A_A 



d t \a = 



-2A+(A_ + (N c + N f )A VA ) + N f A_A ff 
+AvaAct + ^A CT 2 |, 

(10) 



2A ff -4w 4 /i F 1 B) ' 4 [6C 2 (N c ) ff 2 A 



6ff 2 A+] 



1 ,(FB),4 
"4*^1,2 



24 - 9N^ 



-8 U4 4 F) ' 4 {2N c A2-2A_A CT -2N f A CT AvA-6A + A CT }, 



d t A V a = 2A V a - ^v A lfi ] A 



J-.9 2 Ava - 3# 2 A- 



(11) 



(FB),4 



24-3N 



c „4 



N c 



-8 V4 4 F) ' 4 {-(N C + N f )A 2 , A + 4A_A VA - |N f A 2 }. 

Here C 2 (N C ) = (N 2 - 1)/(2N C ) is a Casimir operator of 
the gauge group, and V4 = l/(327r 2 ). For better readabil- 
ity, we have written all gauge-coupling-depcndcnt terms 
in square brackets, whereas fcrmionic self-interactions are 
grouped inside braces. The threshold functions / charac- 
terize the regulator dependence [l4| . For the "optimized" 
linear regulator ^5j, these read 



(F),4 



1 ,(FB),4 _ , _ VA 

2' ^ ~ L 6 : 



;(FB),4 
'1,2 



VA 

6 ' 



(12) 



Next, let us turn to the running of the gauge coupling. 
Even though a non-perturbative estimate can already be 
obtained within the current simple truncation |l(ll |. a re- 
liable estimate to higher loops requires the inclusion of 
many further vertex operators in the truncation [5] In- 
stead, we construct the /3 function according to the fol- 
lowing line of argument: First, we note that gauge invari- 
ance in the RG scheme can be monitored with the aid of 
regulator-dependent Ward-Takahashi identities 0, 0] . 
For the present system, these identities constrain the 
four-fermion contributions to the running gauge coupling 
to be of the form [Tlj : 



f]g2g 2 .r]A, f]-[j) and r\\ are the gluon, quark and quark- 
gluon- vertex anomalous dimensions, respectively, defined 
by Vi — ~c?ilnZi; 77^, in fact, vanishes in the present 
truncation using the Landau gauge. The " symbol de- 
notes those standard contributions which arc not gener- 
ated by the modified Ward-Takahashi identities; for in- 
stance, the difference between 771 and rji is exactly given 
by the term ~ dt\ in Eq. I|13f) . 

Instead of computing f\ g i in the RG scheme, we use 
the four-loop result obtained in the MS scheme 0, 0] 
as an estimate. For a qualitative discussion, let us list 
the first two coefficients (all four known coefficients are 
used in the numerics), 



11 

A) = yN c 



fN f> ft 



34N 3 + 3N f - 13N 2 N f 
3^ ' 



(15) 



34N 

with /3i reversing its sign for Nf > 13N :sl 3 - Since we 

do not expect that the RG- and MS-scheme j3 functions 
exactly coincide, this four-loop approximation introduces 
an error that is estimated below. 3 It should be stressed at 
this point that some crucial properties of the f3 function 
are not scheme dependent; in particular, the existence 
of zeros, i.e., fixed points, is a universal result which is 
relevant for the present problem. 

One further technical point needs to be mentioned: 
whereas the running gauge coupling is dominated by rj g 2 , 
the threshold functions I of the fermionic sector require 
the separate knowledge of tja- As a rough but sufficient 
estimate, we simply employ the one-loop relation for the 
Landau gauge, f? A loop = ^fj g ~2°° p ■ Since the anoma- 
lous dimensions turn out to remain numerically small of 
0(0.1) and the threshold functions only depend on 
this additional higher-order approximation has little ef- 
fect on the result. 4 



III. RESULTS 

The onset of %SB is detected by the four-fermion 
part of our model (2). Chiral symmetry remains pre- 
served in the pointlike limit, as long as the fermion self- 
interactions remain finite. In turn, x^B is signaled by a 
diverging A,. This seeming Landau-pole behavior points 



dtg = fig 2 9 ~ 4u 4 /i 



(F),4 



9 



1 -2w 4 4 F) ' 4 Ec J A 



c a = 1 + N f , c + = 0, c_ = -2, c V a 
where 



-2N f , 



V =Va- 2r?i + 2?^ 

= -4v i9 2 (A) + Pi(2v i9 2 ) + (3 2 (2v ig 2 ) 2 



(13) 



(14) 



is the standard gauge-coupling anomalous dimension, 
trivially related to the standard (3 function by f3 2 = 



3 Note that the additional terms ~ dt\i in Eq. 1131 already con- 
tribute to the resulting (3 function at two-loop order. By approx- 
imating fj g 2 by the MS result, we introduce a potential double 
counting of diagrams at this and higher order. However, the dt\i 
contributions vanish at the fixed point. In addition, wc would 
like to stress that we are using a mass-dependent regularization 
scheme, and universality of the second coefficient in the j3 func- 
tion does not hold and further contributions at this order are 
expected. 

4 As a check, we have used the pure one-loop result rfc oop , which 
is non- vanishing at the fixed point. The corresponding results lie 
well within the estimated error (shaded region in Fig. 151. 
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FIG. 1: Sketch of a typical j3 function for the fermionic self- 
interactions \i\ at zero gauge coupling, a = (solid curve), 
the Gaufiian fixed point Xi — is IR attractive (the second 
fixed point at Ai > corresponds to the IR repulsive criti- 
cal coupling of NJL type). For small a > (dashed curve), 
the fixed-point positions are shifted on the order of a 2 . For 
gauge couplings larger than the critical coupling a > a cr (dot- 
dashed curve), no fixed points remain and the self-interactions 
quickly grow large, signaling x^B. We emphasize that both 
fixed-points values remain finite until the fixed points even- 
tually vanish at a cr . 

to the formation of chiral condensates, since the A^'s 
are inversely proportional to the mass parameter of a 
Ginzburg-Landau potential for the order parameter in a 
(partially) bosonized language, Xi ~ 

Therefore, let us first consider the Xi flow separately. 
For vanishing gauge coupling, the flow is solved by van- 
ishing Aj's, which is the Gaufiian fixed point. This 
fixed point is IR attractive, implying that these self- 
interactions are RG irrelevant for sufficiently small bare 
couplings, as they should be. Since the A^ flows are 
quadratic in the A^'s corresponding to a parabola, there 
are 15 further non-Gaufiian fixed points (2(# Ai=4 ' = 16 
in total) . These are related to the set of critical couplings 
Aj jCr beyond which the system is in the universality class 
of' NJL type, see Fig. □ 

For fixed small gauge coupling, the Gaufiian fixed point 
is shifted a bit to fixed point values of the order A*^ ~ g 4 , 
effectively describing scattering of masslcss quarks. This 
fixed point is still IR attractive, see Fig ^ hence there 
is no xSB as long as the gauge coupling remains small 
enough. 

If the gauge coupling becomes larger than a critical 
value a cr = g 2 T /4ir, all fixed points of the Aj flows are 
removed, such that the Aj's quickly run into a divergence, 
signaling X SB fifr. 

In Fig. |3 we study the four-fcrmion system at fixed 
gauge coupling. For various numbers of colors N c , we 
plot the critical gauge coupling a cr = g^/Air at which 
the fermion system becomes critical and exhibits xSB. 
For an increasing number of colors, a cl - decreases. The 
dependence on the number of flavors is rather weak. 

Now, let us switch on the running of the gauge cou- 
pling. The above-mentioned fixed points in the fermionic 
system become quasi-fixed points which persist for suffi- 
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FIG. 2: Critical coupling for the four-fermion system. From 
top to bottom, the number of colors increases from N c = 2 to 
N c = 8 (thick line corresponds to N c = 3). The red/dark-grey 
(green/ light- grey) line shows the fixed-point gauge coupling 
for N c = 3 at four (two) loop. At the crossing, the critical 
number of flavors can be read off. 

ciently weak gauge coupling; their values arc modulated 
by the logarithmic increase of the gauge coupling. If the 
gauge coupling grows larger than the critical value, the 
fermionic quasi-fixed points vanish (the parabola drops 
below the Xi axis in Fig.^), and the system runs into the 
xSB regime. In this way, the IR physics is solely deter- 
mined by the gauge coupling, as expected. Whether or 
not the system ends up in the xSB regime finally depends 
on the maximal value of the gauge coupling. For a large 
number of flavors, Nf < Nf ' f ', this value is given by the 
IR fixed-point value a*, induced by the higher /3 func- 
tion coefficients. Lowering Nf increases a*, eventually 
exceeding a CI required for xSB. 

We determine the resulting critical number of flavors 
N" by comparing the fixed point of the gauge coupling, 
a* = a*(Nf), to a cr . The task of solving the complete 
coupled set of fixed-point equations dtg 2 = 0, dtX = be- 
comes simplified by the following observation: Eq. i|13|) 
reveals that the only contributions to dtg 2 involving A 
are proportional <9 t A. By definition this term vanishes at 
any fixed point in the subsystem of four-fermion interac- 
tions. Therefore, the gauge coupling at any fixed point of 
the complete system is determined by fj g 2 = 0, which de- 
pends only on g. In Fig. [2 the fixed-point gauge coupling 
a*(Nf) is plotted for N c = 3 with a two- (grccn/light- 
grey) and a four- loop (red/dark-grey) result for f] g 2. 

The critical number of flavors N| r can now be read off 
from the intersection between this line and the critical 
coupling a cr of the four-fermion subsystem for N c = 3 
(thick solid line). In the same way, we determine the 
critical number of flavors as a function of N c , as plotted 
in Fig. 5 



5 An equivalent way to determine Nf r (N c ) is to directly solve the 
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The regulator dependence of the flow equation offers 
a possibility for estimating the error of our quantitative 
results which is introduced by the truncation. Exact re- 
sults for physical quantities like N" are universal with- 
out any regulator dependence. However, this universality 
does not hold for the truncated flow, such that the reg- 
ulator dependence of would-be universal quantities di- 
rectly translates into an error of the truncation. Let us 
start with the fermionic subsystem, for which the point- 
like approximation of vertices and the use of "classical" 
propagators represents the most severe approximation. 
If these truncated nontrivial momentum dependencies of 
vertices and propagators were important for the determi- 
nation of the universal Nf , we would expect a strong reg- 
ulator dependence within the point-like approximation. 
This is because the flow equation is localized in momen- 
tum space and different regulators probe the vertices and 
propagators at different momentum shells; consequently, 
erroneously truncated momentum dependencies of the 
vertices generically imply strong variations of would-be 
universal quantities for different regulators. In this sense, 
the study of regulator dependencies is an important if not 
crucial test of our point-like approximation in the fermion 
sector for the determination of Nf. In the present case, 
the regulator dependence occurs in the threshold func- 
tions I. Our pointlikc truncation is reminiscent of a 
derivative expansion, for which a flow-optimization pro- 
cedure has led to the construction of a linear regulator 
fl5l us ed in Eq. (|I 2(1 . Following the stability arguments 
of [la], we consider the results from this regulator as 
our "optimized" predictions. The opposite limit of non- 
optimized regulators for pointlikc interactions is marked 

by the sharp cutoff [2l| with 4 F) ' 4 = 4 F i B) ' 4 = l ^2 )A = 1, 
neglecting the anomalous dimension. Therefore, a con- 
servative error estimate is given by the difference between 
the sharp cutoff and the optimized regulator. 6 This re- 
sults in the green/dark-grey shaded region in Fig. 
The resulting error is actually very small; in particular, 
we observe strong nontrivial cancellations between the 
contributions of different threshold functions to the final 
result. Since this exactly matches with the expectation 
for universal quantities, it provides evidence for the relia- 
bility of the pointlike truncation for the present purpose. 
As a further evidence for the consistency of the point- 
like approximation in the symmetric phase, we note that 
the gluon as well as the quark anomalous dimension van- 
ish at the IR fixed point in our approximation, implying 
that the use of "classical" propagators is self-consistent 



flow equations lsl- 1131 with Aj(t = 0) = and a small value 
g 2 (t = 0) as initial conditions. One can then find N"' by looking 
for the infimum of all Nf for which the four-fermion couplings do 
not diverge at a finite t, (or alternatively, for the supremum of 
those which do diverge). 

The sharp cutoff can be considered as the limit of a sequence 
of powerlike regulators, the results of which all lie within the 
estimated error. 
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FIG. 3: Critical number of flavors for SU(N C ) gauge theory. 
The result based on the four-loop MS beta function is denoted 
by black circles which lie almost on top of stars, representing 
the three- loop result; black boxes correspond to the two-loop 
beta function. The (inner) green/dark-grey shaded region 
around the four-loop result displays our error estimate for 
the fermionic sector. The (outer) turquois/light-grey shaded 
region shows the approximate gluonic error, estimated by a 
variation of the higher- loop coefficients. The red/grey line 
shows the number of flavors above which asymptotic freedom 
is lost. 



here. Of course, we should stress that the point-like ap- 
proximation is expected to break down for many other 
questions, in particular, those related to the properties 
of the xSB phase. 

The largest uncertainty arises from the gauge sector. 
As a error estimate, one might use the difference between 
the various loop orders in the /3-function contribution r\ g i . 
As depicted in Fig. the difference between two and 
three loops in the MS scheme is of the order AN" ~ 1. 
The difference between the three- and the four-loop re- 
sult is already much smaller. This suggests that, for our 
purposes, the gauge sector is reliably approximated by 
the four-loop /3 function. However, as discussed above, 
the (known) MS result is not expected to coincide exactly 
with the (unknown) result to be obtained with the lin- 
ear regulator used so far. For a quantitative estimate of 
the resulting error, we can perform a comparison at two- 
loop level, where the mass-dependent character of the 
flow-equation regularization already introduces scheme 
dependencies. On a quantitative level, this introduces 
differences on the 10% level For the 3- and 4-loop 

coefficients we assume larger uncertainties on the 30% 
and 50% level, respectively. This results in an error de- 
picted by the turquois/light-grey shaded region in Fig. 
13 Quadratically adding both uncertainties, the resulting 
error is dominated by the gluonic sector. 



7 For an exact computation of the unive rsal two-loop coefficient 
within the functional RG approach, see \ > i . 
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It has been argued |2J| that a regularization scheme 
exists in which the two-loop (5 function for the gauge 
coupling is exact. We stress that this does not imply that 
the two-loop curve in Fig.^is exact, since this particular 
unknown regularization scheme may involve a regulator 
with strong quantitative influence on the fermionic sector 
(possibly requiring operators beyond the pointlikc limit). 
Hence, it is well possible that the resulting curve of that 
special scheme is again close to the estimated four-loop 
result in our scheme. 

It is worth pointing out that the total estimated er- 
ror is significantly smaller than the distance to the line 
Nf" '(N c ) where asymptotic freedom is lost. This strongly 
confirms the existence of a phase of many-flavor QCD 
which is asymptotically free, but has no xSB and no strict 
confinement in the IR. 



IV. SUMMARY AND CONCLUSIONS 

We have used flow equations to determine the criti- 
cal number of flavors N| r separating the phases with and 
without xSB. Our truncation includes all pointlike four- 
fermion interactions allowed by gauge and chiral sym- 
metry which are generated by ladder as well as non- 
ladder processes; thereby, our result represents a sig- 
nificant improvement over conventional approximation 
techniques. Gauge symmetry is monitored by modi- 
fied Ward-Takahashi identities, resulting in an additional 
four-fermion contribution to the flow of the gauge cou- 
pling. We have demonstrated that these terms have no 
effect on N| r , owing to the structure of the modified 
Ward-Takahashi identity (by contrast, these terms are 
of quantitative importance, for instance, for the Landau- 
pole problem of QED H3). 

Our findings confirm the existence of a sizable pa- 
rameter region N" < Nf < Nf' f ' where chiral sym- 
metry remains unbroken, in agreement with earlier re- 
sults 0, IS El El- Our result for the critical num- 
ber of flavors, for instance, in SU(3) gauge theory, is 
N^ r = 10.0±0.29(fermion)+J-^(gluon). The errors re- 
sult from the fermionic and gluonic sectors of our trun- 
cation, respectively, the quantitative influence of which 



we have estimated from the regulator dependence of uni- 
versal quantities. 

Our restriction to a truncation in terms of quark and 
gluon fields with pointlike interaction does not facilitate a 
study of the x>3B regime; for this, the inclusion of bosonic 
composite fields along the lines of [l3l I20L 12(1 W\ repre- 
sents a powerful tool. 

Nevertheless, important aspects of the unusual nature 
of the phase transition can already be read off from 
the fermionic sector: as long as the system is in the 
symmetric phase, N" < Nf < Nf ' f ', the fermionic self- 
interactions approach fixed-point values in the IR which 
are finite numbers. Even arbitrarily close to the phase 
transition, these numbers remain finite and arc elements 
of a compact region in parameter space, Aijfe^o € By 
contrast, the self-interactions diverge in the broken phase 
where all fixed points have vanished. As a consequence, 
an infinite region in parameter space outside Q\ remains 
inaccessible. Varying Nf continuously across the phase 
transition coming from the symmetric side, the fermionic 
self-interactions jump discontinuously at Nf = Nf r . Now, 
the Aj's are inversely proportional to the mass parameter 
of a Ginzburg-Landau potential for the order parameter 
in a (partially) bosonized language, Xi ~ 1/mf; hence, 
our observation in the fermionic sector agrees with the 
observations of - IS El that the phase transition is not 
conventionally second order, even though the chiral or- 
der parameter is known to change continuously across the 
chiral phase transition |l3l |. Our work thus reveals fur- 
ther aspects of the nature of this type of zero-temperature 
phase transitions which may find further application in 
related systems such as QED3 or models of clcctroweak 
symmetry breaking. 
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